Dynamical decoupling and noise spectroscopy 
with a superconducting flux qubit 

Jonas Bylander^, Simon Gustavsson^, 
Fei Yan^, Fumiki Yoshihara'^, Khalil Harrabi^^ George Fitch^, 
David G. Cory^'^'^, Yasunobu Nakamura'^'^, Jaw-Shen Tsai'^'^, William D. Oliver^'^ 

^Research Laboratory of Electronics, ^Department of Nuclear Science and Engineering, 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA. 
^The Institute of Physical and Chemical Research (RIKEN), Wako, Saitama 351-0198, Japan. 
''MIT Lincoln Laboratory, 244 Wood Street, Lexington, MA 02420, USA. 
^Institute for Quantum Computing and Dept. of Chemistry, University of Waterloo, ON, N2L 3G1, Canada. 
^Perimeter Institute for Theoretical Physics, Waterloo, ON, N2J, 2W9, Canada. 
^Green Innovation Research Laboratories, NEC Corporation, Tsukuba, Ibaraki 305-8501, Japan. 
^^Present address: Physics Dept., King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia. 

January 26, 2011 



Abstract 

The characterization and mitigation of decoherence in natural and artificial two-level systems (qubits) is funda- 
mental to quantum information science and its applications. Decoherence of a quantum superposition state arises 
from the interaction between the constituent system and the uncontrolled degrees of freedom in its environment. 
Within the standard Bloch-Redfield picture of two-level system dynamics, qubit decoherence is characterized by 
two rates: a longitudinal relaxation rate Fi due to the exchange of energy with the environment, and a trans- 
verse relaxation rate r2 — ri/2 + T^p which contains the pure dephasing rate T^p. Irreversible energy relaxation 
can only be mitigated by reducing the amount of environmental noise, reducing the qubit's internal sensitivity to 
that noise, or through multi-qubit encoding and error correction protocols (which already presume ultra-low error 
rates). In contrast, dephasing is in principle reversible and can be refocused dynamically through the application 
of coherent control pulse methods |[T}{3|. In this work we demonstrate how dynamical-decoupling techniques can 
moderate the dephasing effects of low-frequency noise on a superconducting qubit |4-6| with energy-relaxation 
time Ti = l/Fi — 12 /is. Using the CPMG sequence [7j^ with up to 200 7r-pulses, we demonstrate a 50-fold 
improvement in the transverse relaxation time T2 over its baseline value. We observe relaxation-limited times 
_ 23 /is w 2 Ti resulting from CPMG-mediated Gaussian pure-dephasing times in apparent excess 
of 100 IIS. We leverage the filtering property of this sequence in conjunction with Rabi and energy relaxation 
measurements to facilitate the spectroscopy and reconstruction of the environmental noise power spectral density 

(PSD) imiioi. 
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Several multi-pulse sequences developed within the field of nuclear magnetic resonance ^ (NMR) have re- 
cently been applied to mitigate noise in qubits based on atomic ensembles fTTI, semiconductor quantum dots 
1.1 2 , 1 3 1 , and diamond nitrogen-vacancy centres [.14.15il . We extend these methods to the realm of superconducting 
quantum devices, and subject a remarkably long-lived qubit to varying levels of longitudinal and transverse noise 
by rotating the qubit's quantization axis, against which we characterize the baseline coherence rates Fi, and 
r^. We evaluate three dynamical-decoupling pulse protocols: the Carr-Purcell |7| (CP); Carr-Purcell-Meiboom- 
Gill Is] (CPMG); and Uhrig dynamical-decoupling lfT6l (UDD) sequences. The narrow-band filtering property of 
the CPMG sequence enables us to sample environmental noise over a broad frequency range 0.2-20 MHz, and we 
observe a l//"-type spectrum which we independently confirm with a Rabi-spectroscopy approach. We further- 
more characterize the environmental noise from 5.4 to 21 GHz by monitoring the qubit's relaxation rate ll9l [T0]| . 

The device is a persistent-current qubit (Figs.[T^ and Al), an aluminium loop interrupted by four Al-AlO^j-Al 
Josephson junctions. When an external magnetic flux $ threading the loop is close to half a superconducting flux 
quantum $o/2, the diabatic states correspond to clockwise (counterclockwise) persistent currents Ip = 0.18 /iA 
with energies zLhe/2 = ±/p$b, tunable by the flux bias <i>b = $ — ^o/'^- At <i>b = 0, the degenerate persistent- 
current states hybridize with a strength hA ~ h x 5.3662 GHz (Figs.[TJ) and A2b), where h, = h/2Tr and h is 
Planck's constant. The corresponding two-level Hamiltonian is HIS] 

n^-^[{e + de)a, + {A + 6A)a,], (1) 

which includes noise fluctuation terms 6e and SA, and a^^z are Pauli operators (Fig. A2). The ground (|0)) and 
excited (|1)) states have frequency splitting ajgi = Ve^ + A^ and are well isolated owing to the qubit's large 
anharmonicity, uji2/^qi ~ 5. The environmental noise leading to fluctuations 6e (e.g., flux noise) and 6 A (e.g., 
critical current and charge noise) physically couples to the qubit in the e - A frame (equation [T]i. However, their 
manifestation as longitudinal noise (dephasing) or transverse noise (energy relaxation) is tunable |17| by the flux 
bias $b and determined, respectively, by their projections Suj^' onto the qubit's quantization axis a^' (which makes 
an angle 9 = arctan(£ / A) with a^) and 6uj±' onto the plane perpendicular to (t^/ . 

The chip is mounted in a ^He/^He dilution refrigerator with 12-mK base temperature. For each experimental 
trial, we initialize the qubit by waiting sufficient time (^ 1 ms) for it to relax to its ground state. We drive the 
desired quantum-state rotations of angle Q by applying calibrated in-phase (Xq) and quadrature (Ye) harmonic 
flux pulses to the qubit loop. The pulses comprise Gaussian envelopes with a typical standard deviation cr = 1.2 ns 
and truncated at ±3<t. The qubit readout has 79% visibility (Fig. A3) and is performed in the energy basis by 
determining the switching probability Psw of a hysteretic dc SQUID, averaging over several thousand trials (see 
Appendix). 

We begin with a spectroscopic characterization of our device and its baseline coherence times. The qubit 
level splitting ljqi is measured via saturated frequency spectroscopy (Fig. [TJ)), and at low-power it exhibits a 
Lorentzian full-width-at-half-maximum (FWHM) hnewidth A/(fwhm) = 0.18 MHz at $b = (Fig.[TJ;). The 
energy relaxation is generally exponential and its time constant Ti — 12 ± 1 /is (Fig. [TJl) is remarkably long among 
superconducting qubits |6|, a feature we leverage in this work in conjunction with quantization axis tunability. We 
observe similarly long decay times at $b = for the Hahn spin-echo, T^^^e = 23 /is (Fig. [TJl), Ramsey free 
induction, T2 — 2.5 /xs (Fig.[TJ), and Rabi oscillations, Tr = 13 /is (Fig.[TJ). Although the spin echo and Rabi 
exhibit an apparently exponential decay at $b = (they are essentially Ti -limited at this flux bias), in general, their 
decay functions are non-exponential. Furthermore, the dephasing times decrease (rates increase) dramatically away 
from $b = (Figs.|2}; and[3j)) due to the qubit's increased sensitivity to the dominant (5e-noise (flux noise) in this 
system fTSl (increased \du!oi/de\), which, as we will demonstrate, can be mitigated with multi-pulse dynamical- 
decoupling sequences to push coherence times up towards the Ti -limit. 

Decoherence in superconducting qubits has been studied theoretically (]9]|T7][l9]|20l and experimentally ifTOl 
[17][T8 |. Each noise source A is characterized by its PSD, which quantifies the frequency distribution of the noise 
power, S\{uj) = {l/2n) (A(O)A(t)) exp(— iwi). In contrast to a simple Bloch-Redfield picture, which 

presumes weakly coupled noise sources with short correlation times Tc <C Ti, T2 resulting in purely exponential 



2 



a bed 




Flux detuning, <t|^ {mftg) Frequency, f (GHz) t (|.is) 





^ Ramsey: ^ 


h ft i 1 L I % t Rabi oscillations Tp^lS^s ^^'^ 


1 








r2"=2.5|.is 


1 1 1 1 



2 4 6 5 10 15 20 25 

^ (l^s) X (|JS) 

Figure 1 : Qubit device and characterization, a. Device and biasing schematic: An aluminium super- 
conducting loop interrupted by Josephson junctions (crosses) with a read-out dc SQUID, b, Frequency 
spectroscopy of the qubit's |0) — > |1) transition, c. Spectroscopy at $b = (arrow in b). d. Echo 
decay (blue triangles) and relaxation from the excited state (black dots) at <i>b = 0. In the insets, t is a 
time delay and Xq symbolizes a rotation of the Bloch vector by the angle O around the axis a-^- The 
red squares indicate the time of read out. e. Free-induction decay (Ramsey fringe) at $b = 0. f, Rabi 
oscillations at <l>b = 0. 



decay functions, more generally the decay functions are non-exponential for noise sources with long correlation 
times, singular near w « (e.g., 1 //-type noise at low frequencies, relevant to dephasing). 

A superposition state's accumulated phase (p{t) ~ {ojQi)t + S(p{t) diffuses due to adiabatic fluctuations of the 
transition frequency, S(p{t) — (duiQi/dX) dt'6X{t'). For noise generated by a large number of fluctuators that 
are weakly coupled to the qubit, its statistics are Gaussian. Ensemble averaging over all realizations of the stochas- 
tic process SX{t), and taking the sources A to be independent, the dephasing is {exp[i 6(p{t)]) = exp[~XN{t)], 
with the coherence integral 

XN{r) =r^Y. (^) ' ^^(^) 9n{^^ 

where t is the free evolution time and will denote the number of vr pulses in the pulse sequences lfm[T6ll . 
Equation dU expresses that the PSD S\{ljj) is filtered by a dimensionless weighting function giq determined by 
the pulse sequence, and the aggregated A-noise translates to dephasing through dujQi/dX, the qubit's longitudinal 
sensitivity to A-noise. For Gaussian noise with spectral distribution S\{uj) = A\/uj at low frequencies, the 
coherence integral results in a Gaussian decay function, xn{t) = (Fi^r)^. 

The Ramsey free-induction and Hahn spin-echo dependence on flux bias (Fig. |2]:) are both apparently con- 
sistent with Gaussian distributed, l//-type noise. Ramsey free induction, the free evolution of a superposition 
state for a time t (Fig. with no tt pulses), has a filter function go peaked at w = (Fig. |2j3) and is sensitive 
to low-frequency longitudinal noise Sll!z'{uj — > 0). Inhomogeneities in the precession frequency luqi from one 
reaUzation of the pulse sequence to the next lead to a decay of the averaged signal. We denote such fluctuations 

"quasi-static" noise and characterize them by a noise variance af^2 /^^a" duj S\{ljj), with cut-off frequencies 

W;^. and cjyy determined respectively by the averaging time over all trials and the typical free-evolution time during 
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a single trial. In contrast, the Hahn spin-echo sequence |[T], a single tt pulse applied at time t/2 (Fig. with one tt 
pulse), has a filter function gi peaked away from a; = (Fig.|2}3) and is less sensitive to quasi-static noise. We plot 
the decay rates 1/Te versus flux bias $b for Ramsey and echo in Fig.|2}; (see also Fig. A4), where, for purposes of 
comparison amongst different decay envelopes, Tg parameterizes the time T2 to decay by a factor 1/e independent 
of the exact decay function (see Appendix). At 4>b = 0, 5A-noise is the dominant longitudinal noise that limits 
the Ramsey decay, yet it is refocused with a single 7r-pulse resulting in the Ti -limited exponential echo-decay in 
Fig.[TJl. As |$b| is increased, both the Ramsey and echo decay rates increase due to the qubit's increased longi- 
tudinal sensitivity to fc-noise (see Figs.[2t and A4). The (5e-noise is too large for the echo to refocus efficiently, 
due to its high-frequency tail, and we find best-fit phase decay functions that are Gaussian, x(r) — (F^ f(e) r)^. 
We extract the ratio F(^'b)/r,^,E($b) ~ 4.5, as expected for 1// noise ifTTl . with an equivalent flux-noise 
amplitude flS'l A$ = (1.7 /i$o)^- Importantly, we note that in this analysis and related works 11711181 . the PSD 
was presumed a priori to take the form 1//" with a = 1. 

Numerical simulations (see Appendix), including the measured Ti decay at each $b, linearly coupled and 
uncorrected quasi-static noises, and uncorrected dynamic 1// noise from 10^ to 10^° Hz, reproduce the entire 
$b-dependence of and T2 e using the parameters in Table 1 and are also consistent with equation (|2]). 

In order to mitigate higher-frequency noise than what the Hahn spin-echo can efficiently refocus, we further 
shape the filter function by applying additional tt -pulses (Figs. and A5). The filter (7iv(w, r) depends on the 
number N and distribution of 7r-pulses lil lul6..21.,22J during the total sequence length r. 



where Sj £ [0, 1] is the normalized position of the centre of the jth 7r-pulse between the two 7r/2-pulses, and Ttt is 
the length of each tt -pulse. As the number of pulses increases for fixed r, the filter function's peak shifts to higher 
frequencies (Fig.|2]3), leading to a reduction in the net integrated noise (equation m for l//"-type noise spectra 
with a > 0. Alternatively, for a fixed time separation r' = t/N (valid for N > 1), the filter asymptotically peaks 
near uj' = 27r/4r' as more pulses are added. In principle, one can adapt the filter function to suit a particular noise 
spectrum via the choice of dynamical-decoupling protocol. 

We have evaluated three different dynamical-decoupling protocols relevant for 1 //-type power law noise spec- 
tra ||23l - |25l . The CP and CPMG sequences |7, 8 1 are multi -pulse extensions of the Hahn echo with equally spaced 
TT-pulses whose phases differ from the initial tt/2 pulse by 0° (X^) and 90° (Y^), respectively (Fig.|2^). The UDD 
sequence lfT6l has -pulse positions defined by Sj = sin^ ( 2N+2 ) ■ 



In Fig. 12};, we include the 1/Te decay rates for CPMG dynamical-decoupling sequences with TV = 2 . . . 48 tt- 
pulses along with the Ramsey (N = ) and Hahn echo (N ~ 1) already discussed. The decay rates monotonically 
improve towards the l/2ri-limit as the number of 7r-pulses increases, extending the range around $b = for 
which dephasing is negligible. Outside this range, as the qubit's sensitivity to Se (flux) noise grows with increasing 
l^bl, increasing the number of 7r-pulses mitigates the noise sufficiently well to achieve a desired decay rate. 

At a specific flux bias $b = 0.4 m$o^ where the qubit is highly sensitive to 6e noise, the CPMG sequence 
gives a marked improvement in the decay time up to w 200 7r-pulses (Fig.|2}l), beyond which pulse errors 
begin to limit the CPMG efficiency. We achieve a 50-fold enhancement of jj^CPMC ^ygj. jj^g Ramsey T2 (Fig.|2}l), 
and well over 100-fold improvement in the Gaussian pure dephasing time T^p. The CPMG sequence performs 
about 5% better than UDD, indicating that the 1// (5£-noise spectrum exhibits a relatively soft (if any) ultraviolet 
cutoff l,21.i25 1. and it dramatically outperforms CP, as -pulse errors appear only to fourth order with CPMG, 
whereas with CP, Xt^ errors accumulate to second order ll26l . 

We use the filtering property of the CPMG sequence to characterize the (5e-noise spectrum. The filter (7Ar(w, r) 
is sufficiently narrow about uj' that we can treat the noise as constant within its bandwidth B and approximate 
equation Q as xn{t) ~ (dujoi/d^)"^ Se{uj') gwi^^' ,t') B (Fig. A6). We compute uj' and B numerically 
for each TV and t used in the CPMG measurements of Fig. |2]l. The measured decay function contains three 
decay rates: dephasing F^^'' and exponential relaxation ri/2 during the total free-evolution time r, and pulse- 



5Ar(w,r) = - — — 1 + (-1)-^+^ exp(za;r) + 2^(-l)^ exp(ia;5jr) cos(a;r^/2) 



(3) 
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induced decay Fp during Nt-j^ (we assume Fp to be independent of A^). Conceptually, Figs. |2^-g illustrate for 
= 84 how we divide out the Fi and Fp (assuming exponential pulse-induced decay) components from the raw 
data (e) to obtain the Gaussian phase decay (f), and then compute ^^(a;) (g). More rigorously, we only use that 
method to determine a starting point at a single frequency, and then use a recursive method to obtain the remainder 
of the spectrum without presuming a functional form for the pulse-induced decay, the dephasing, or the noise 




Figure 2: Dynamical decoupling pulse sequences, a-b, CPMG pulse sequence and filter 
function gN{uj,T) for a fixed total pulse-sequence length r. c. Decay rates (inverse of 1/e 
times) vs. flux detuning: Free induction (green squares) and CPMG (coloured dots) with N = 
1, 2, 4, 6, 8, 10, 16, 20, 24, 30, 36, 42, 48. SoHd lines are calculations using equation ^ with param- 
eters in Table 1. d-g. Measurements at $b — — 0.4m$o (£/27r — 430 MHz), d, 1/e decay time 
under A^-pulse CPMG, CP, and UDD sequences. The simulation (red line) assumes perfect pulses, 
e. Population decay under an 84-pulse CPMG sequence. Gaussian fit with fixed exponential decay 
contributions Ti = 12 /is, Tp = 1.75 /is. f. Phase decay of the signal in e after dividing out Ti and Tp. 
g, Noise PSD calculated from the data in f. 



5 



spectrum (see Appendix). Both approaches yield nearly identical l//"-type spectra with a slight increase in the 
measured PSD above 2 MHz; we plot the recursively extracted PSD S^iuj) over the region 0.2-20 MHz in Fig.|4l 
Interestingly, the PSD estimated in this manner is better approximated by a 1//" power law |27 1 with a — 0.9 < 1 
(solid, red line) with noise amplitude — (0.8 /i$o)^^ obtained by fitting the lower-frequency, linear portion of 
the PSD. Projecting this line to higher frequencies comes within a factor 2 of the transverse noise at frequency A 
as extracted from the relaxation measurements described below. 

We confirmed the spectrum over a similar frequency range by analyzing the decoherence during driven evolu- 
tion, which provides an independent means to rotate the quantization axis with respect to the noise sources when 
viewed in the rotating frame |17,28|. A transverse driving field at frequency uj results in Rabi oscillations with 
angular frequency flu = y^fl'^ + (Auj)'^ « il + (Aaj)^/2f2, where Aw = w — cjoi. Integrating the oscillations 

over a normal distribution with variance a^, we obtain the quasi-static decay function C(t) = (l + (ut)^) 

where u = (e/wpi)^ o'e/^- Along with Fi, the noise at the Rabi frequency r[^' — nSx{QB.) comprises the usual 
exponential Rabi-decay rate 

r„.(|r, + ir<f')co.,'.+ (-l-)"lr<-, (4, 

where cos^ 61 « 1 as the quantization angle is small. The combined decay function is ({t) exp (— Frt) (Fig. [3^; 
see also Appendix). 

To determine S'A(r2R), we measured the Rabi oscillations vs. $b with fixed Rabi frequency ilR. For each 
l^R, we find the £-independent part of the rate by fitting the envelope of the oscillations at e = 0. The rate 
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Figure 3: Decoherence during driven dynamics, a, Rabi oscillations with f2R/27r ~ 2 MHz at 
£/2tt = 225 MHz. The red line envelope is a fitting using (^{t) and Fr. The black line shows the Fr 
decay only, and the green line the (^(i) -envelope contribution, b, Rabi-decay rate Fr vs. flux detuning 
at r2R/27r = 2 MHz, with a paraboHc fit to equation (HI used to obtain Si.{fl^) — ^ r[^^ 
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was too small to distinguish accurately from Fi, consistent with its correspondingly small quasi-static noise cr^ 
(Table 1). Then, for e 7^ 0, we divide out the known quasi-static contribution ({t) and fit to the parabolic term in 
equation (|4|i, from which we obtain T^^'^ (P^g-Wp)- Using this approach, we find S^{fl) to be consistent with the 
1//" noise obtained from the CPMG measurements (Fig.]?]). 

We now turn to transverse noise, i.e. (5e-noise at e = 0, and (5A-noise at e ^ A (Fig. |4] inset), at the qubit 
frequency woi responsible for energy relaxation Fi. In the low-temperature limit, /cbT <C hujoi, where the 
environment cannot excite the qubit, the golden-rule expression for Fi in a weakly damped quantum two-level 
system is 

ri = f E (^ys.M^^sM, (5) 

with duj±' /dX the qubit's sensitivity to transverse noise and S{luqi) the total PSD (see Appendix). 

We apply a long (^ Ti) microwave pulse to saturate the transition and monitor the energy decay to the ground 
state, using equation ^ to determine S{ujoi ) over the frequency range A < u>oi < 27r x 21 GHz by tuning $b- At 
$b = (cjoi = A) and using a measurement-repetition period t^op > 1 ms, we observe Ti = 12 ± 1 /is as shown 
in FigllJl. As t^cp becomes shorter than 1 ms, the decay becomes increasingly non-exponential, which we attribute 
to the residual presence of non-equilibrium quasiparticles generated by the switching SQUID during readout. We 
observe structure in the Fi data due to environmental modes (e.g. cavity modes, impedance resonances) with 
uncontrolled couplings to the qubit (Fig.|4]i. For comparison, we plot the expected Johnson-Nyquist flux noise in 
Fig.|4]due to the R= 1/G = 50 fl environment mutually coupled with strength M = 0.02 pH to the qubit via the 
microwave line, 

('")=2^(,a$j 1-e-W^-^- 

This known noise source falls about 100 times below the measured PSD at /oi ~ A/2tt (red dot in Fig.lH where 
the relaxation is due solely to 6e noise, indicating that the dominant source of energy relaxation lies elsewhere. 
The crossover /c between the effective 1//- and /-type flux noises (Fig. |4| occurs between A/27r and k-gT/h, 
where T — 50 mK is the approximate electronic temperature of our device 11011301 . 



Table 1 : Quasi-static noise parameters used in simulations, and coherence times. 
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Ti 


T2 
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$b = m(E>o 


12 ^s 


2.5 /is 


23 ^ls (N = 1) 


9 


$b = 0.4 m(f>o 


12 ^s 


0.27 us 


13 us (N = 200) 


48 



In Ramsey-fringe and Hahn-echo simulations, we describe the Gaussian noise distributions by their 
standard deviations, a\, obtained by integrating the 1// noises over the bandwidth given by the ex- 
perimental protocol (cut-off frequencies uj-^ and a;^^, see text). At e = 0, the dephasing improvement 
under a Hahn echo is greater than the theory would suggest for (5A 1// noise that extends to high 
frequencies; the lower uo^ gives consistency (see Appendix). The equivalent flux and normalized 
critical-current noise amplitudes. Ax, are values derived from the Ramsey and echo data assuming a 
power law 1//" with a = 1 and that all noise in e and A is flux and critical-current noise, respec- 
tively; they are consistent with previously reported values 1*181 1291 . Using these parameters in simula- 
tions yielded agreement with A^-pulse dynamical-decoupling data, consistent with equation (|2]l. The 
coherence times are given at two bias points dominated by 5/S. and Se noise, respectively. 
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The dynamical-decoupling protocols demonstrated in this work comprise the same types of simple pulses that 
are used for quantum gate operations and therefore require little additional overhead to implement. Integrating 
refocusing pulses into qubit control sequences, e.g., by forming composite gates that incorporate both quantum 
operations and refocusing pulses, will lead to lower net error rates in systems limited by dephasing 1311 . Despite 
observing levels of 1// flux and critical-current noise similar to those observed ubiquitously in superconducting 
qubits and SQUlDs [i6ji29J, we could mitigate this noise dynamically to increase the pure dephasing times beyond 




Frequency, f (Hz) 



Figure 4: Noise-power spectral density (PSD). Multi-coloured dots, ^e-noise PSD (0.2-20 MHz) 
derived from CPMG data at <i>b — — 0.4m<i>o (see text). Colours correspond to the various N in 
Fig Eh-c; grey dots for data up to iV 250. Yellow squares, Je-noise PSD (2-20 MHz) derived 
from Rabi spectroscopy (Fig.[3l see text). Diagonal, dashed lines. Estimated 1// flux (red) and SA 
(blue) noise inferred from the Ramsey and echo measurements (cf. Ax and in Table 1). SoUd, red 
line. Power-law dependence, Se{f) = A^/(27r/)", extrapolated beyond the qubit's frequency, A/27r. 
The parameters A'^ = (0.8 /i$o)^ and a = 0.9 were determined by fitting the low-frequency, linear 
portion of the CPMG PSD data before the slight upturn beyond 2 MHz (see Appendix). The shaded 
area covers a ± 0.05. Green dots. High-frequency e and A PSD inferred from energy-relaxation 
measurements above A/2tt — 5.4 GHz. Purple line. Guide to indicate linearly increasing Nyquist 
(quantum) noise, including the eigenbasis rotation (see inset); dots indicate transverse Se (red) and 
SA (blue) noises. Inset, Graphic representation of the quantization axis (grey arrows of fixed length) 
with the qubit's (Z') eigenstate tilted from the "laboratory" frame (Z) by the angle 9. Fields e(fE'b) 
and A point in the X and Z directions, respectively. Red and blue double-arrows indicate transverse 
noise. 
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0. 1 ms, more than a factor 10^ longer than the intrinsic pulse length. However, despite having a remarkably long 
energy relaxation time Ti = 12 fis, the transverse relaxation T2 ~ 2Ti was ultimately limited by it. Dynamical 
decoupling protocols go a long way to refocusing existing levels of 1// noise and achieve long coherence times, 
and the main emphasis is now on identifying and mitigating the noise source(s) that cause energy relaxation. We 
note for further study that the PSD power law obtained experimentally by the CPMG technique, when extended to 
higher frequencies, falls within a factor two of the measured transverse noise, suggesting the possibility that the 
microscopic mechanism responsible for low-frequency dephasing may also play a role in high-frequency relax- 
ation. 



Appendix 
Measurement set-up. 

We performed our experiments at MIT, in a dilution refrigerator with a base temperature of 12 mK. The device was 
magnetically shielded with 4 Cryoperm-10 cylinders and a superconducting enclosure. All electrical leads were 
attenuated and/or filtered to minimize noise. 

We used the Agilent E8267D microwave source, and employed the Tektronix AWG 5014 arbitrary waveform 
generator to create 1/Q modulated pulses, and to shape the read-out pulse. 



Description of the qubit. 

We fabricated our device at NEC, using the standard Dolan angle-evaporation deposition process of Al-AlO^j-Al 
on a Si02/Si wafer (Fig.lABi). 

The persistent-current, or flux qubit ||4l|5][32l[33l consists of a superconducting loop with diameter d ^ 2 fim, 
interrupted by four Josephson junctions (Fig. lAll) . Three of the junctions each have the Josephson energy — 
210 GHz, and charging energy Ec = 4 GHz; the forth is smaller by a factor a = 0.54. The ratio of energy scales 
puts the device in the flux limit, E'j / Eq ~ 50, thus making the phases across the Josephson junctions well defined. 
The geometric and kinetic loop inductances are negligible compared to the Josephson inductance: 

Lg-^orf--2pH, Li, = ^loXll/S ^ SOpU, Lj = $o/27r/c ~ lOnH, 

where we used Al = 100 nm, I — 10 /im, and S" = 20 x 250 nm^. 



a 








1 |im 




dc + (j,w 



Figure Al: a. Scanning electron micrograph (SEM) of a device (qubit and SQUID shown) with iden- 
tical design parameters as the one measured during this work 1 18 1. b. Schematic. Qubit loop (shaded) 
and galvanically coupled read-out SQUID. The crosses are Josephson junctions; R, bias resistors; C, 
shunt capacitances; L, inductances. 
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When the external magnetic flux $ threading the qubit loop is close to half a magnetic-flux quantum, $o/2, 
the qubit's potential energy exhibits a double-well profile with quantized energy levels ±?i£/2 — ±/p$b, where 
$b — $ — $0/2 is the flux bias and Ip — 0.18 /iA the persistent current. The ground states in the left and right 
potential wells correspond to the diabatic states \L) and \R) of opposite circulating persistent currents. A tunnel 
coupling TiA/2 between the right- and left- well qubit states opens an energy gap, hA, between the ground and 
excited states, |0) and at flux degeneracy, $b — 0. The two-level Hamiltonian - analogous to a spin- 1/2 
particle in a magnetic field - is in the "laboratory" frame 

n = -^[{e + de)a, + {A + 6A)a,], (7) 

where hujoi = 'hsj + is the energy-level splitting; Oj are the Pauli matrices; and be and are the noise 
fluctuations. The quantization axis makes an angle Q = arctan(e/A) with (and an angle it — 6 with the 
persistent-current eigenstates), so that the qubit is first-order insensitive to flux noise when biased at the "sweet 
spot" £ = 0. (We write the Hamiltonian in this way rather than with swapped x and z indices as in several 
previous papers, so that our pulses will be along X and Y in agreement with the habitual NMR language.) It can, 
alternatively, be written in the qubit's eigenbasis, 

n = -^h {ujoi^z' + ScJz'CTz' + Soj±>d-±') , (8) 

see Fig. IA2fa . Here d-±i denotes that the transverse spin component can include both ax' and ay/ . The longitudinal 
and transverse noises, SlUz' and Suj±/, are further described below. 



Z = Z' 




Y=Y' 



X = X' 




{X,Y,Z}\s the laboratory frame 
{X,' Y',Z'] is the qubit's eigenframe 



> 

T3 
"ft 
X 
(J 

o 
-c 




0.48 



0.50 



0.52 



Figure A2: Two-level Hamiltonian in different eigenbases, and simulated energy spectrum. 

a. The Hamiltonian (|7]i in the "laboratory" frame {X, Y, Z} has the e field (flux-bias) along X, per- 
pendicular to the plane of the qubit loop, and the fixed tunnel coupling A along Z. The Hamiltonian 
(O in the qubit's eigenframe {X', Y', Z'} makes an angle 9 with Z. The figure shows both frames 
under two bias conditions (e = [6* = 0°] and e = 4A), with the laboratory frame fixed in space. The 
two frames coincide when e = 0. b. Simulated energy spectrum. 
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Description of the read-out SQUID. 

Our hysteretic dc-SQUID ||34l . see Fig. lAll has critical current 1^ — 4.5 /iA; normal resistance Rjy = 0.25 kf7; 
mutual qubit-SQUlD inductance Afq-s — 21pH; on-chip capacitors C ^ lOpF, inductors L ^ 0.1 nH, and 
bias resistors Rj = 0.2 k51 and Ry — 1 k57; and further cold and room-temperature resistors and filters. A 
coil with mutual coupling M^c — 0.6 pH to the qubit provides the dc flux bias, and an on-chip antenna with an 
estimated M^w ~ 0.1 pH the microwave excitation. The shunt capacitors bring the plasma frequency down to 
Wp/27r = 2.1GHz. 

The qubit's magnetization, resulting from the persistent currents, modulates the SQUlD's switching current. 
For qubit read out, we apply a sample-and-hold current pulse to the SQUID (Fig. IA3fa ). and use a threshold 
detector, after room-temperature amplification, to register the presence or absence of a voltage, conditioned on the 
qubit being in state \L). We determine the switching probability, Pgw, statistically by repeating this measurement 
several thousand times. 

The pulse is produced by a digital arbitrary-waveform generator (Tektronix AWG5014) with 250 MHz analog 
output bandwidth. Due to enhanced relaxation when a current is flowing through the SQUID, a rapid sample pulse 
is important for good read-out visibility. A hold current enables the room-temperature electronics to register the 
voltage pulse before retrapping occurs, but a low and short hold pulse limits the quasiparticle generation, which is 
also beneficial for the visibility. We obtain an optimal read-out visibility of 79% (Fig. lA3b ). 

In addition to biasing the SQUID, the read-out pulse couples flux into the qubit with the effect of shifting 
the states adiabatically, |0) \L) and |1) -> \R), before the measurement takes place, so that the states can be 
distinguished by the projective measurement in the persistent-current basis. 

The optimal SQUID dc-bias current /b with respect to noise coupled into the qubit was very close to /b = 0, 
indicating highly symmetric SQUID junctions. In this set of experiments, we therefore did not apply any dc bias 
to the SQUID while manipulating the qubit, cf. Ref. ITSl . 




1.54 1.58 1.62 1.66 
SQUID sample voltage, (V) 



Figure A3: a. Read-out pulse (voltage across a l-kTl bias resistance). The optimized pulse shape 
is programmed to have a 4-ns sample pulse, including 1-ns rise, 1-ns sample, and 2-ns fall times, 
followed by a 3-/is hold plateau at 17.5% of the sample voltage, b. Read-out visibility at <l>b = 0. 
Scans of the SQUID's switching probability, Psw, vs. the height of the sampling pulse. We obtain 
79 % read-out visibility between the qubit's ground and excited states. Relaxation during read out 
leads to an imperfect determination of the excited state (17 % dark counts at the optimal Vs). 
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Sensitivity to noise. 

Following the approach of Ithier et al. {VT\, we will evaluate the noise terms in the qubit's Hamiltonian, 

i-L — i^oi'^z' + 6llIz'(Jz' + Suj^'a^i) [equation (|8])]. 

We consider both flux and critical-current noise (affecting A), and express the noises as sensitivity derivatives that 
translate the noise fluctuations to a change in the Hamiltonian, i.e. 5\ — ^ SijJ. At $b — 0, the first-order noise 
(5A is much larger than second-order 5e (flux) noise, as confirmed through simulation. It is therefore sufficient to 
expand to first order to explain our Ramsey and echo data, 

5uj^, — 5\ + . . . and Suj^i ~ 6X + . . . . (9) 
o\ oX 

Energy relaxation is related to noise that is transverse to the qubit's quantization axis, Suj±r (A), at the frequency 
of the level splitting wqi- 

Pure dephasing, on the other hand, is related to low-frequency fluctuations of the qubit's energy-level splitting. 
To evaluate the longitudinal first-order term in equation ([9]) we use the chain rule. 



dujQi dujQi de (9woi d/S. 

By geometry. 



(10) 



oe loqi o/\ wqi 

From spectroscopy measurements we infer the e sensitivity to flux noise, A = $, 

de 

e=^ = 27rx l.lGHz/m$o, (12) 

while A is insensitive to $ noise. 

Since he = 2Ip $b and /p = Ic\/^ — 1/(2q;)^, e is sensitive also to critical-current noise, A — Iq, 



de e 



(13) 



For the A sensitivity to fluctuations in = 5Ic/Ic, a numerical simulation gives dA/dIc — 2?: x 8 ■ 10 Hz/A, 
and with Ic = 0.4 /iA we obtain 

dA 

Ki = — = 27r X 3.01 GHz. (14) 

die 

Taken together, the first-order fluctuations are therefore 

^SX^ ^^5^, + ^^51. + A,,. (15) 
OA cjQi woi ic ^ai 

Flux noise dominates, except very near e = 0, where the SA noise is needed to account for the observed decay. 
(The second term is negligible.) 
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Decoherence. Ramsey and echo phase decays — 1/ /-flux noise. 



In the Bloch-Redfield formalism, valid for weakly coupled, short-correlated noise, the dynamics of two-level 
systems is described by the longitudinal and transverse relaxation rates, Fi ~ 1/Ti and r2 = I/T2, respec- 
tively [ 17 28 |. We assume that the qubit is coupled to many fluctuators, which, in concert and independent of their 
individual statistics, will yield a Gaussian noise distribution due to the central limit theorem. In systems where the 
noise is regular at the frequency of the qubit's energy-level splitting, luqi, relaxation and dephasing factorize. The 
longitudinal relaxation is exponential due to the many uncorrected contributions of transversally coupled noise 
at ujQi. The pure-dephasing rate, T^p, associated with low-frequency, quasi-static, longitudinally coupled noise 
(inhomogeneous broadening), combines with Fi to give the rate r2 = Ti/2 + T^. However, the exact form of 
the time dependence of the dephasing component is determined by the noise-PSD. (That is, when the dephasing is 
non-exponential, the inverse time constant can strictly no longer be interpreted as a rate.) 

At $b = 0, the echo decay is nearly Ti limited in our device, and therefore practically indistinguishable from 
an exponential. Relaxation contributes by Fi = 1/2 Ti = 43 x lO'^/s or 6.8 kHz to the low-power spectroscopic 
line width A/(fwhm) ~ 0.18 MHz. This exceeds the expected A/(fwhm) = l/'''^2 ^y only 0.05 MHz, indica- 
ting little power broadening, given the free-induction decay rate, Fp = I/T2 ' with the measured = 2.5 /^s. 

Biased away from (f>b, we find a Gaussian spin-echo phase decay, exp[— (F;p e t)^], consistent with a 1// flux- 
noise model ifTTI (singular PSD at a; = 0). As the Ramsey fringe decays considerably faster it is hard to distinguish 
between exponential and Gaussian decays. The bias dependencies for echo and Ramsey decays are 

r<^,E/F(^'b) = (^<I.??E/f)^ 

where the numerical factors 77 differ due to the different echo and Ramsey filtering functions: rj^p^p = \n{l/Ldi,-t) 



dujoi 



(16) 




Flux detuning, O^, (mOg) Flux detuning, (mOp) 

Figure A4: Decay rates vs. flux detuning, a. Energy -relaxation rate, Fi, and Gaussian phase-decay 
rates for the Hahn echo, Fi^ e, and Ramsey free induction, T^^f, after subtracting the exponential Fi 
decay. Straight, black lines are fits to F;^ E/F(*i'b) for flux noise only; see the text and equation ( fT6l l. 
b. Blow up of the data in the dashed-box region in a, along with simulated Ramsey- and echo- decay 
rates including both 6e and SA noises. For the echo, the grey line is the rate obtained for a 1// noise 
in 6 A with the 1-MHz ultra-violet cut off given by the experimental protocol. The red line is the 
rate obtained with the lower ultra-violet cut off: uj^^/2tt < 0.1 MHz. The exact rate depends on the 
detailed shape of the noise cut off, which is unknown. 
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and ?7<p.E = In 2. Their ratio, r<p.F('i'b) /ri^.E(*f b) ~ 4.5, is in accordance with our data, and the magnitude. 



of the noise S{uj) = agrees with previous results [18J, see Fig. lA4l 

Numerical simulations of Ramsey and echo. 

The (E>b dependencies of T2.E are reproduced in numerical simulations. We simulated the Ramsey fringe 

by numerically solving the Bloch equations, including the measured Ti decay and linearly coupled quasi-static 
noises, averaging over many realizations. Each run had stochastic, normally distributed Se and 5 A deviations from 
their average values: The 1// flux noise, ^(a;) = gives a^j^-K = lOMHz when integrated from 1 Hz to 

1 MHz; noise in A dominates near = 0, and we get agreement for cta/Stt = 0.06 MHz, uncorrected with the 

flux noise and obtained from 5'a('^) = ^a/I^^I with Aa = (4 • 10^^)". 

In the echo simulation we additionally took into account dynamic noises from 10^ to 10^° Hz. We obtained 
the noisy time series for e and A by inverse-Fourier transforming the amplitudes of the 1// noises with random 
phases for each Fourier component, and then evaluated the Schrodinger evolution operator in discrete time. 

Near <I>b = 0, the experimental data shows a greater echo-improvement in T2 than the simulation would suggest 
for these parameters (see Fig. lA4b ). This discrepancy would be explained by a lower ultra-violet A-noise cut off, 
with faster decay than 1// above ~ 0.1 MHz. 

Numerical evaluations of the coherence integral, equation ( [20b below, agree with our simulations. 

Pulse calibration. 

We calibrate the rotations to < 1 % accuracy by applying a tune-up sequence of pulses, akin to methods used in 
NMR [35,1 . A rigorous measurement of gate errors should be done with randomized benchmarking, for example. 

Dynamical-decoupling pulse sequences. 

Collin et al. |[36ll and Ithier et al. ifTTl have employed some NMR methods beyond the Hahn spin echo to ma- 
nipulate a superconducting qubit. In this work, we apply multi-pulse, dynamical-decoupling pulse sequences to 
significantly enhance the coherence times and to facilitate spectroscopy of the environmental noise. 
During the Carr-Purcell-Meiboom-Gill (CPMG) multi-echo sequence Il7][8l, defined as 



and illustrated by the rotations in the Bloch sphere (Fig. lASI ). the transverse component of the Bloch vector is refo- 
cused along the axis of the refocusing pulses (here Y), whereas the perpendicular component (X) is randomized. 
CPMG is a development on the Carr-Purcell (CP) sequence Q, which has identical pulse positions, but where all 
rotations are along the same axis. The third pulse sequence that we investigated, the Uhrig dynamical-decoupling 
(UDD) sequence l,16.,37J . has normalized pulse positions defined by 



(17) 




(18) 




(19) 



The dephasing under a certain sequence is described by the coherence integral. 




(20) 
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The dimensionless filter function gN{^, t) depends on the number and distribution of tt pulses 111 Il|16ll2ni22l[37]| . 



(c.r)2 



(21) 



i + (-i) 



N 

i+^exp(iwr) + 2^(- 



-1)^ exp{iLd5jT) cos(a;r^/2) 



(22) 



where Sj G [0, 1] is the normalized position of the centre of the jth tt pulse between the two 7r/2 pulses. The 
TTT-dependent factor assumes square pulses, but is a reasonable approximation in our case. Alternatively to equa- 
tion {TH . one can define a function Fn{uj,t) — |yAr(tj,T)p that filters the phase noise, S{uj)/lu^, as in, e.g., 
Refs. |[IT][l6l|2T]|22l|321. 

The Ramsey free induction (N — 0) and the Hahn spin echo (N — 1) have weighting functions goiuj, t) — 
sinc^(ajT/2) and gi{u!, r) = sinc^(a;T/4) sin^(cjT/4), respectively. Note that .gi(0) = 0, indicating a suppression 
of the low-frequency part of the noise for the Hahn echo. 

The CPMG sequence is inherently robust to field inhomogeneities when the following criteria are fulfilled ||26l : 
(i) The effective rotation axis is oriented in the transverse (XY) plane (true when the driving frequency is resonant 
with the level splitting); (ii) the magnetization (Bloch vector) is initially aligned with the rotation axis (Y); and 
(iii) the rotation angle is tt. Errors that occur due to deviations from (i-iii) can be quantified by computing the 
propagator over one cycle of the sequence. For an initial state Y, one finds that F^r errors appear only to fourth 
order (CPMG), whereas errors accumulate to second order (CP). 




•^^11/2 



-►X 




UDD sequence 

_U_J \ L 



N TT-pulses 

J L 



W TT-pulses 




Sequence repeated N times 

CPMG sequence 

■^71/2 '^K ^JI ^71 ^7t ^71 ^JI ^JI ^jr ^JI ^ir ^7l/2 

I I I I I I I I I I I I- i 

W jT-pulses 

CP sequence 

Xk/2 ^^^'"^^^^^ ^ ^ ^k/I 

_LJ \ \ \ \ \ \ \ \ LU 



J \ LJU 



Figure A5: Dynamical-decoupling pulse sequences. Rotations of the Bloch vector during the CPMG 
sequence (pulses along X and Y in the laboratory frame.) Timing of the CPMG, CP, and UDD pulse 
sequences for iV = 10. 
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Determination of the PSD from the CPMG-decay data. 

The qubit is subjected to different decay mechanisms during the CPMG pulse sequence. We identify three decay 
rates, and fit the qubit's population to the function 

Psw{r) = Pa + a exp {-Fit) exp(-rp(Tp)) exp (-^^(t)) , (23) 

keeping the constants a and Pq fixed for all N. 

• Relaxation. For the duration of the sequence, there is energy relaxation with the constant rate Fi = 1/2 Ti, 
where Ti = 12 /is. 

• Decay during pulses. In addition, there is pulse-induced decoherence during the total pulse time Tp = 
Ntt^ + 2 T3^/2- (For experimental reasons, instead of 7r/2 pulses we often use the equivalent 3tt/2 pulses.) 
We do not need to assume that this decay has any particular form, e.g. exponential, as it will be divided 
out of the calculation. However, a simple exponential decay with Fp = 1/(1.75 /is) gives a good overall fit 
independent of N. Although much higher than the Rabi-decay rates described in the manuscript, this rate is 
reasonable given that we observe an increased Rabi-decay time at high driving amplitudes. 

• Dephasing. The fit parameter F^^-* is the dephasing rate during the total free-evolution time r under an 
A^-pulse CPMG sequence, which provides us with a means of figuring the noise PSD, S'(aj). Only the PSD 
appears in the decay function, because the noise statistics were taken to be Gaussian, and all information is 
included in the second-order correlation function. 



Dividing equation (|23l l with itself for two free-evolution times, ri and T2, we get 



q{Ti,T2) = ^ =exp(-Fi[ri -r2]) cxp ( - [xn{ti) ~ xn{t2)]) ■ (24) 

Psw[T2) - Pq \ /V y 



For a pulse separation r' = t/N (and large N), the filter (|2TI) is peaked near 



- §J- (25) 



It is narrow enough, that we can treat \u)'{t'), t' j as a delta function peaked at w', and regard the noise in the 
qubit's transition frequency as constant within its bandwidth. Aw, as illustrated in Fig. |A6b . This provides us with 
a means of characterizing the noise at a frequency u by varying r. We can rewrite the integral in equation ( |20l ) as 

J^duj gj\j(^uj'{T'),T'^ Auj Jp°°d w — w'), and do the approximation 

XN{r) ^ (^^y r's(u;\r')) gr,(c.'ir'),r')Ac.. (26) 



We know Auj{t), and (a;'(T'), r' j from a numerical calculation for each N. 

We can use equation ( |26l ) to evaluate the second factor on the right-hand side of equation ( l24b . 

- '^l S ((^' iTl)j gN,l[i^'{T2),T2^ AuJl - T2 S (u' {T2)j ^AT (^w'(r2), T2^ Au2 

(27) 

This gives us the difference between the noises at uj'{ti) and Ld'{T2), so that we can find one if we know the other. 
To get the absolute noise level, we have to assume S{uj) for some value of uj' = uj'{t'), e.g., a value consistent 
with the PSD found in the spin-echo measurements, equation ( fTTI i. or indeed choose xn = 1 for r = and go 
from there. Following this procedure, we obtain a number of points representing the noise PSD over 0.2-20 MHz. 
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The points within the range 0.2-2 MHz lie tight together and are the least susceptible to error, which is why we fit 
this data to a straight line, rendering the reported 1//" dependence: A'^ — (0.8 /i$o)^ and a ~ 0.9. Integrating 
the noise PSD from these parameters yields = 8 MHz (compared to 10 MHz) for the quasi-static noise, slightly 
changing the slope of the calculated F^ p in Fig. IA4I which presumed a = 1. However, the slope of the PSD is 
not necessarily constant over many orders of magnitude, which could explain the slight deviation for the Ramsey, 
being sensitive to low-frequency noise. 

Alternatively, we can fix the rates Fp, Fi, and then fit the remaining, Gaussian decay xn{t) — (T^t)'^. Then 
each measured point, Psv,{t), is mapped onto S{uj'{t)) by using equation ( l26l l. However, the using a Gaussian 
decay function is equivalent to an a priori assumption of a = 1. 

Figure lA6b shows a numerical evaluation of the coherence integral (l20l i. assuming a noise spectral density 

S{uj) = A/uj. (28) 

The calculation shows, that this method underestimates the decay by about 5 %, since it disregards the part of the 
noise that falls within the harmonics of g^, the first of which occurs at 3a;' and which is about 10 dB lower than 
the main lobe at oj'. 




2n/4x 10 20 30 40 50 

Angular frequency, co Total free-evolution time, t (|is) 



Figure A6: a. The filter function is narrow, and the noise PSD is smooth enough, that we can obtain the 
noise PSD at the filter's centre frequency by evaluating the coherence integral XNi^), equation ( |26] |. 
b. Fitted phase decay, exp[— XAr(T)], under an iV = 84 7r-pulse CPMG sequence with ~ 6.8 ns, at 
the flux detuning e/27r = -430 MHz ($o = -0.4m3>o). Blue is the function exp{-{r^^K)'^) with 
a fitted F^^^ = 8.5 /is. Red is the integral ( |20] | with S{uj) — A/lu. Note the excellent agreement 
between the blue and red lines. Dashed is the expression ( |26] i with the same parameters as for the red 
line. 



Decoherence during driven evolution. 

We obtain information about the flux noise in the 2-20 MHz frequency range also by analyzing the decoherence 
during driven evolution. 

Driven by a near-resonant, transverse field Hi = —Ti il{t) cos(cji) a^, the static part of the Hamiltonian in the 
rotating frame is 

n^^^hiAojaz + na^), (29) 
where Aw = u — wqi- As a result, the qubit oscillates at the Rabi frequency 

fiYi = y^n^ + (Aw)2 + {Auj)^/2n (30) 

around the effective field, which makes an angle r/ = arctan(Aa;/f2) with the qubit's (T^ axis. 
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The decay of the Rabi oscillations has contributions from noise at the frequencies ujqi and 17 IfTTl . and from 
the low-frequency, quasi-static flux noise with variance cr^. 

Static noise. We first describe the decay due to quasi-static noise. Flux deviations from the bias point e lead to a 
frequency detuning Aoj = {dujoi/de)Ae = (£/a;oi)A£, normally distributed with a standard deviation a^, 

7V(A£) =(771/2 exp(-(Ae)V2a2). (31) 

Integrating over N{Ae), the Rabi oscillations decay as 

= J d{Ae) N{Ae) cos {fin t + (/)). (32) 

Under the approximation Auj ^ fl, the Fresnel-type integral (equation[32]i becomes IfTTl 

C(t) = (1 + (utf)^^^'^ cos (^^t + ^ arctan(Mt)j 

where 



(33) 



o \22 / \22 



de J Q \ujQi J il 



(34) 



Noise at the Rabi frequency. We now turn to the exponential decay of the oscillation envelope (r2 in the rotating 
frame), 

Tr^ ^Ficos^^+irn, (35) 

where the component that depends on the noise at the Rabi frequency il^ has two contributions which we denote 

r[?'andr(-': 

Tn = TrS^iQK) sin^ + nSAi^R) cos^ ^ = ^n^ sin^ 6* + ^It^ cos^ 0. (36) 
We can approximate cos^ « 1, for e <C A, and write equation (|35] l as 



^n = -^i + 7;^r +{ — ] 7;^'n'- (37) 



Fitting procedure. In order to determine S{il^), we measured the Rabi oscillations vs. "I>b, with il^ fixed 

1/4 

for each set of data. The combined decay envelope from Eqs (|33] |37] | becomes (l + (wt)^) x exp(— rpi). 

We observe that, at e = 0, the model predicts an exponential decay with the rate |ri + For a given 

JIr, we find this rate by fitting, and then keep it fixed while fitting the Rabi envelopes vs. e. At low amplitude 
(r2R/27r ~ 2MHz), Fr nearly reaches its limit (3/4)ri. Dividing out the known quasi-static contribution, and 
then fitting the envelope to the parabolic equation (|37] |. we obtain the flux-noise dependent from which we 
can calculate S'e(il). 



Energy relaxation; high-frequency spectroscopy. 

Figure lAT] shows the variation over repeated measurements of Ti relaxation. To explain the high-frequency noise, 
leading to Ti relaxation, we model the e and A channels of the environment as ohmic resistors. The microwave an- 
tenna's designed mutual inductance to the qubit is M — 0.02 pH with i? = 50 termination. The Johnson-Nyquist 
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Figure A7: Ti relaxation. Repeated measurements of Ti relaxation from the excited state, at flux 
degeneracy = 0, showing how Ti varies by about 10% from one measurement to the next. The 
size of the error bars from the exponential fits are on the order of the variation. Each trace took 1 min 
40 s to measure; Ti may vary on a shorter time scale than that. 

(thermal and quantum) noise of an environment at T 50 mK (a reasonable effective electronic temperature on 
the chip) i^ 

Sfi-)--(?-yM^^^^, (38) 
^ ^ ^ 27r \d<P J 1 _e-W=BT' 

which simplifies to 

5JN(^) ^ ^ (^|£^ {2k^T + 2n^)/R. (39) 

This expected noise falls about 100 times below the measured PSD at /qi — A/2tt, indicating that intrinsic, 
microscopic noise mechanisms dominate the flux noise. With R unchanged, we would infer an effectively ~ 10 
times larger M to account for the observed noise. 

Within this work, we did not identify which noise sources were responsible for the A noise, which can arise 
from, e.g., critical-current noise or charge noise due to various microscopic mechanisms. 

Taking the tilted quantization axis into account, so that the transverse noise is constituted by e noise at e = 0, 
by A noise at e A, and by combined projections of the two noises in-between, we arrive at the combined 
expression for the effective transverse Johnson-Nyquist noise, 

5a (0.01) = ^5^(0.01) + "^"S ^' ^a'(^oi). (40) 

As mentioned above, this effective noise is obviously greater than the expected extrinsic noise sources (Johnson- 
Nyquist as well as 1// noise). The subgap conductance, typically several hundred to several thousand times the 
normal resistance, may limit the achievable Ti -relaxation times lf38l . a topic worth further studies. 
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